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Comprehensive catecholaminergic projectome 
analysis reveals single-neuron integration of zebrafish 
ascending and descending dopaminergic systems 
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Essential components of animal behaviour are modulated by dopaminergic (DA) and 
noradrenergic circuitry. In this study, we reveal at cellular resolution the complete set of 
projections ('projectome') of every single type of DA and noradrenergic neurons in the central 
nervous system of zebrafish larvae. The most extensive DA projections are established by 
posterior tubercular ofp-dependent neurons, with individual somata integrating the ascending 
DA system, the descending diencephalospinal, as well as the endohypothalamic circuitry. 
These findings suggest a major role in the modulation of physiology and behaviour for 
ofp-dependent DA neurons, which correlate with the mammalian A11 group. We further 
identified an endogenous subpallial DA system that not only provides most of the local 
DA projections, but also connects to the ventral diencephalon. The catecholaminergic 
projectome map provides a framework to understand the evolution and function of these 
neuromodulatory systems. 
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Distinct dopaminergic (DA) neuronal groups establish local 
and far-ranging projections to modulate circuits controlling 
many essential aspects of behaviour, cognition and physi- 
ology. Much attention has focused on the DA systems, given their 
clinical relevance in neurological and psychiatric disorders such as 
Parkinsons disease. Knowledge of the complete set of DA projec- 
tions would improve the understanding of neural circuit function in 
the normal brain, the pharmacological treatment of diseased states 
and the requirements for neurorestorative therapies of degeneration. 
However, the complexity in the number of neurons and connections 
in the vertebrate brain has so far prohibited global approaches to 
comprehend connectomes'. 

Studies in rodent, bird, amphibian and fish vertebrate models 
have revealed conserved and evolutionary dynamic aspects of cat- 
echolaminergic (CA) systems 1,2 , which in the central nervous system 
(CNS) consist of DA and noradrenergic (NA) neurons 1 . An advan- 
tage of zebrafish embryos and larvae as an experimental system is 
their translucent CNS, small enough to be imaged at cellular reso- 
lution in toto. Furthermore, a functional nervous system is estab- 
lished after only 4-5 days of embryonic development and enables 
zebrafish larvae to perform complex behaviours such as swimming 
and hunting 3 . Most of the analysis of zebrafish CA systems is based 
on the detection of Tyrosine hydroxylase (TH) protein or th tran- 
script expression, which starts before the end of the first day of 
development 4 7 . A second paralogous th gene, th2, is expressed later 
than th, and is predominantly observed in caudal hypothalamic DA 
neurons during early larval stages 810 . Zebrafish DA neurons form in 
the diencephalon and telencephalon, but not in the mesencephalon, 
wherein groups A8-A10 form in mammals 4 " 6,810 . However, ventral 
diencephalic DA neurons of the zebrafish posterior tuberculum that 
emanate ascending projections to the subpallium (SP) have been 
suggested to provide functions similar to the mammalian mesodi- 
encephalic DA systems 11,12 . Mammalian DA groups are numbered 
in a caudal (Al) to rostral (A17) order 1,2 , and may be compared with 
zebrafish DA systems on the basis of anatomy and molecular mark- 
ers (see Supplementary Note 1). Similar to All DA cells in mam- 
mals, zebrafish ventral diencephalic DA groups DC2 and DC4-6 are 
specified by the transcription factor Orthopedia (Otp) 1315 . The other 
DA groups correspond to their mammalian anatomical counterparts 
in the hypothalamus (A 12, A14-A15) and ventral thalamus (A 13). 
Similar to mammals, NA cell groups are located in the locus coeru- 
leus (LC), medulla oblongata (MO) and area postrema (AP) 16-18 . 

Although CA axonal tracts in larval 4,6,12,19 and adult 5,11,16,18,20 
zebrafish have been described, projections of individual neurons 
and their target areas have not been identified previously. In this 
study, we investigate the complete projections ('projectome' 21 ) of DA 
and NA systems in the larval zebrafish. To accomplish this goal, we 
tagged individual CA neurons by green fluorescent protein (GFP) 
expression using genetic mosaics and generated a large data set 
from which the projection behaviour of each type of th- expressing 
CA neuron in zebrafish larvae could be deduced. Our data revealed 
that Al 1-type DA neurons are the major far-projecting DA neurons 
in zebrafish, with single neurons projecting to telencephalon, dien- 
cephalon, hindbrain and spinal cord. Furthermore, we showed that 
a local subpallial DA system provides the major DA projections into 
the anatomical equivalent of the striatum in zebrafish. Our findings 
establish an anatomical framework for understanding the DA mod- 
ulation of behaviour in zebrafish. 

Results 

Genetic labelling of single CA neurons and their processes. 

We introduced a Ras membrane-anchored EGFP (rasGFP) 
reporter at the ATG of the reading frame into a bacterial artificial 
chromosome (BAC) vector containing a 154kbp upstream region 
and most of the zebrafish tyrosine hydroxylase (th) gene. We chose 
th because at larval stages it is expressed in nearly all zebrafish 



CA groups, whereas the paralogous th2 gene is expressed only 
in a small number of hypothalamic DA neurons 810 . Expression 
of thrasGFP BAC is specific for th (Supplementary Fig. SI). 
thrasGFP was injected into one-cell-stage embryos and efficiently 
labelled individual CA somata and their axons in genetically mosaic 
larvae. Systematic injections established a large data set of labelled 
neurons and enabled recording of axonal projection patterns for 
each population of CA neurons in zebrafish CNS at 4 days post 
fertilization (d.p.f.). At this stage, most CA neuronal clusters and 
tracts observed in the adult brain are established 4,12 . Diencephalic 
DA neuronal axons reach the spinal cord in as little as 12 h after 
onset of th expression 6,12 . Therefore, most DA axons are expected 
to have reached their larval stage targets at 4 d.p.f. Larvae (h = 408) 
coimmunostained with anti-GFP and anti-TH antibodies were 
analysed by confocal microscopy. Anti-GFP immunoreactive (ir) 
cell somata were identified as members of each of the previously 
defined DA neuronal groups on the basis of coexpression of TH, 
thus revealing that thrasGFP drives GFP expression in all known 
THir CA groups. A total of 310 larvae containing 670 GFP-labelled 
somata with distinct, non-intersecting projections contributed to 
the final data set used in this study (Supplementary Table SI). 

We established an image analysis workflow to obtain a data set 
of single CA neuron projections and to identify and classify their 
projection behaviour (Fig. 1; Supplementary Methods, Supplemen- 
tary Movie 1 and Supplementary Fig. S2). This imaging technique 
did not allow us to distinguish between axons and dendrites. Target 
area anatomical identities were confirmed from histological stains. 
To comprehensively analyse the CA projection network, we first dis- 
sected the characteristics of projections emanating from neurons of 
each CA group to include directionality (based on the body axes) and 
type (for example, arbourizing, branching, midline crossing, local; 
Supplementary Fig. S2). For further classification, major projections 
that are proximal to the cell soma and target the contralateral side 
were indicated as contralateral, whereas distal branches that originate 
from major ipsilateral projections but target contralaterally were doc- 
umented as midline crossing processes. The second step was to visu- 
ally trace all CA projections to their anatomical target areas (Fig. 2). 
Statistical analysis indicates that our data set is large enough to iden- 
tify most ( > 80%) of the projection behaviours and target areas. 

Overview and classification of projections and target areas. 

Figure 2 summarizes projection patterns of all CA groups. Our fur- 
ther analysis structured the CA groups into five classes. First, NA 
systems comprise LC, MO vagal and AP neurons that share the 
main feature of far- ranging projections (Fig. 3), whereas NA neu- 
rons of the MO interfascicular zone predominantly connect locally 
(Supplementary Note 2 and Supplementary Fig. S3). The second 
class includes local as well as intermediate length projections of DA 
neurons: the larval DC1 group, which contains ventral thalamic and 
periventricular posterior tubercular DA neurons, and the pretectal 
DA group (Fig. 4; detailed in Supplementary Note 3). Class 3 DA 
groups project exclusively locally or to adjacent brain territories: 
hypothalamic anterior preoptic, preoptic, DC3 and DC7 groups, as 
well as olfactory and retinal DA groups (Fig. 4; Supplementary Note 
4 and Fig. S4). In class 4, DA connections to distant brain regions 
appear to be predominantly formed by the ofp-dependent dien- 
cephalic DA groups 2, 4, 5 and 6 (Figs 5-7). Finally, in the subpal- 
lium local DA neurons provide the predominant DA arbourization 
as well as descending projections to the hypothalamus (Fig. 8). 

Complexity of projection patterns. Comparison of projection 
patterns shown by single neurons within each group led us to note 
that high individual neuron target variability exists. Taking DC 5 
(n = 1 10) as an example, although the neurons predominantly project 
into the endohypothalamic tract (eht) and their descending axons 
terminate in the hindbrain or spinal cord, we identified different 



2 NATURE COMMUNICATIONS | 2:171 | DOI: 10.1038/ncomms1171 | www.nature.com/naturecommunications 

© 201 1 Macmillan Publishers Limited. All rights reserved. 



NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1171 



ARTICLE 




Figure 1 1 Analysis of individual zebrafish larval catecholaminergic neurons and projections in the context of the intact brain, (a) Scheme of 
dopaminergic (DA; blue) and noradrenergic (NA; red) cell groups and catecholaminergic (CA) tracts (grey) in 4d.p.f. zebrafish larvae 12 . Left: lateral view 
of forebrain and midbrain; right: dorsal view of CNS. Dark blue indicates dorsal-most and light blue ventral-most DA groups. Diencephalic DA groups: 
(DC1) ventral thalamus and periventricular posterior tuberculum; (DC2, DC4) large neurons in posterior tuberculum; (DC5, DC6) medium-sized neurons 
in posterior tuberculum and hypothalamus; (DC3) medial hypothalamus; (DC7) caudal hypothalamus; PO, preoptic region; POa, anterior preoptic region; 
Pr, dorsal pretectum; AC, retinal amacrine cells. Telencephalic DA groups: SP, subpallium; OB, olfactory bulb. NA groups: LC, locus coeruleus; MO, 
medulla oblongata interfascicular zone and vagal area and AP, area postrema. ac, anterior commissure; act, anterior CA tract; eht, endohypothalamic 
tract; mlct, medial longitudinal CA tract; pc, posterior commissure; poc, postoptic commissure; poht, preopticohypothalamic tract; prp, pretectal 
projections; prtep, pretectotectal projections, (b-e) In silico separation of individual neurons (arrowheads in b) and their projections from whole-mount 
confocal image stack raw data, (b) Unprocessed maximum intensity projection (MIP) of a 4d.p.f. larval brain. Blue masks D and L were manually drawn 
on the MIPs using ImageJ to separate soma and projections, based on visual discrimination. Bottom-most panel shows an example for computation of 
3D segmentation mask A/I for a single neuron in the AP, based on the drawn masks D and L. The resulting 3D mask M was applied to the green channel to 
isolate the neuron and its projections, (c) Visualization of segmented AP neuron shown in (b) with projections. Quality control (QC; cyan) was carried out 
in 3D space within the region of interest (between dotted lines) to ensure the masks did not omit any GFP-labelled features belonging to the particular 
neuron, (d, e) Visualization of segmented individual soma and projections of GFP-labelled (d) amacrine cell (arrowheads) and (e) NA neuron in the 
MO interfasicular zone (IZ; arrowheads). See Supplementary Movie 1. Anterior at left, dorsal at top; anti-TH (red), anti-GFP (green); Dor, dorsal view; 
Lat, lateral view. Scale bars, 20um. 



combinations of target choice by processes arising from individual 
somata. These variations may be qualitatively assessed from the 
selected examples in Supplementary Fig. S5. A detailed analysis 
of categories of projection behaviours is found in Supplementary 
Fig. S6. We emphasize this observation to highlight the complexity 
of the CA systems even though further investigation into this aspect 
is not within the scope of this study. To assess projection behaviours 
typical for a specific neuronal group, we conducted Fishers exact 



test to mark projection behaviours typical for most neurons of each 
group (Supplementary Fig. S7). 

Hindbrain NA neurons widely project throughout the CNS. The 

LC, MO vagal area (VA) and AP NA groups innervate a broad range 
of target regions extending from the forebrain to the spinal cord. The 
typical projection topology of individual somata starts with a single 
major axon leaving the soma, generating numerous subsequent 
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Figure 2 | Target areas of catecholaminergic projections. The target areas (columns) of projections of each catecholaminergic (CA) neuronal group 
(rows) are indicated quantitatively by the colour code representing the percentage of total neurons per analysed group that project into the area. The total 
number of neurons analysed is shown in the furthest right column. For CA groups in which individual somata extend projections to multiple target areas, 
the percentages add up to more than 100. Lavender indicates the neuronal group is liquor contacting. Projections pass through the commissures with the 
exception of a few cases in which some processes were noted to terminate there: anterior commissure as target of locus coeruleus and area postrema 
noradrenergic neurons; postoptic commissure as target of anterior preoptic, preoptic, DC1, 2, 4 and 5 dopaminergic neurons. 1' (for local projection) 
denotes the termination of neuronal processes within the region in which the soma is located. DC1 includes populations located in the ventral thalamus 
and periventricular posterior tuberculum. 



branches in up to six different brain territories (Figs 2 and 3; Supple- 
mentary Fig. S2) . For LC NA neurons, this often includes ipsilateral and 
contralateral areas (14 of 35 somata; Fig. 3a; Supplementary Movie 2). 
LC neurons predominantly emanate ascending projections as far 
anterior as the SP and olfactory bulb, as well as to the posterior tuber- 
culum and hypothalamus. Telencephalic targets of the LC that were 
not reflected in Figure 2 include the pallium (n = 3) and the pallial- 
subpallial border area (n = 3). In addition to forming local processes, 
LC NA neurons also send descending projections down the spinal 
cord. Of the LC- derived projections that target the anterior hindbrain 
(Fig. 2), -35% terminated in the cerebellum. The main NA contribu- 
tions to the pretectum (11 of 13 NA projections) and tectum (12 of 
13 NA projections) were observed to originate from the LC (Fig. 2). 

For analysis, we subdivided MO into the interfascicular zone, 
VA and AP 18 . The majority of somata from the VA (Fig. 3b) and 
AP (Fig. 3c; Supplementary Movie 3) were observed to send a far- 
reaching axon and multiple processes into the telencephalon and 
diencephalon (41 of 46 MO somata), in addition to strong local 
projections within the MO (Fig. 2). Approximately 69% of these 
MO-VA neurons formed forebrain projections to the SP, posterior 
tuberculum or hypothalamus. Similarly, most (28 of 33) of the exam- 
ined AP NA neurons were found to emanate long ascending axons 
from an individual soma that target either both the telencephalon 
and different ventral diencephalic regions, or only one of these two 
target areas. They also arbourize into or sparsely innervate the AP 
(Fig. 3c). The remaining AP NA neurons projected only within the 
hindbrain (5 of 33). Furthermore, for < 10% of the cases we noted 
minor NA branched processes from the MO that contribute to 
spinal cord projections (Supplementary Fig. S2). 

4 



Projections of Otp-dependent DA neurons. Within the larval DA 
system, DC2, DC4, DC5 and DC6 are major contributors to the eht 
(201 of 280 DA processes to eht) and DA tracts to the hindbrain and 
spinal cord (237of 243 DA descending projections; Figs 2, 5 and 6). 
Most of these neurons develop complex branching patterns, often 
contributing contralateral projections in the hindbrain (49%; Sup- 
plementary Fig. S2). Minor processes originating from DC2, DC4 
and DC5 target the pretectal and anterior tectal areas (Fig. 2). A high 
proportion of individual neurons from these four clusters project 
locally within the ventral diencephalon (that is, local processes 
contribute to poc, eht and caudal hypothalamus), as well as to the 
hindbrain and/or spinal cord. However, the distinct Otp-dependent 
clusters behave differently with respect to the relative contribution to 
these target areas. Both local and descending projections were iden- 
tified for 78% of all examined DC2 neurons, 90% of DC4, 64% of 
DC5 and 49% of DC6 (Fig. 7; Supplementary Figs S2 and S6). 

These findings prompted us to analyse whether any of the four 
Otp- specified, but anatomically distinct, DA groups exhibit special- 
ized projection behaviours. Statistical clustering analysis revealed 
that, overall, descending projections are typical for Otp-dependent 
neurons (Supplementary Fig. S7). However, the neuronal groups 
differ in the caudal extent of their projections. Significant numbers 
of neurons documented for DC2 (45%), DC4 (57%) and DC5 (29%) 
contribute to the spinal cord, whereas those of DC5 (63%) and DC6 
(86%) preferentially ended within the rhombencephalon. The DC5 
group targeted the LC significantly more often than other DC groups. 
In contrast to other otp- expressing groups, the most anterior group, 
DC2, more frequently projects ventrally into the preoptic region 
and laterally into the lateral diencephalon. This suggests that either 
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Figure 3 | Long-range projections of hindbrain noradrenergic neurons, (a-c) Z projections of confocal stacks in dorsal (Dor) and lateral (Lat) views, 
(a) A locus coeruleus noradrenergic (NA) neuron (arrowheads) projects into multiple regions of the brain including posterior tuberculum (PT), anterior 
preoptic region (POa) and the spinal cord (SC) via a descending branched axon. A projection from the soma ascends and crosses the midline via the 
posterior commissure (pc; arrow in middle bottom panel), and forms an arbour of fine processes within the tectum (T; magnified in rightmost panels). 
Middle and right panels are magnification of boxes in leftmost panels. See Supplementary Movie 2. H, hypothalamus, (b) A medulla oblongata vagal area 
(VA) NA neuron (arrowheads) sends a far-reaching axon into the tel- and diencephalon. The long axon forms branches within the posterior tuberculum 
to innervate the hypothalamus as well as the contralateral side via the anterior catecholaminergic tract (act) and the anterior commissure (ac). A short, 
branched process (top right panel) appears to target the rhombencephalon. Local fibres emanate from the proximal part of the axon. Middle and right 
panels are magnification of boxes in left panels, (c) A single area postrema (AP) NA soma (arrowheads) sends long ascending axons terminating in 
contra- and ipsilateral sides of the diencephalon. The right projection ends in the dorsal thalamus (THd), whereas the left innervates the preoptic region 
(PO; magnified in below dor and lat panels). See Supplementary Movie 3. Anterior at left, dorsal at top; 4d.p.f. larvae; anti-TH (red), anti-GFP (green). 
Scale bars, 20 um (d) Schematic overview of NA projections in 4d.p.f. zebrafish larvae. Projection pathways (indicated by lines) and target areas 
(indicated by arrowheads) are depicted in lateral (left) and dorsal (right) views. Short arrows (see MO and AP) denote formation of local arbours or 
processes. For visual clarity in dorsal views, lines and arrows are not depicted left-right reciprocally. Our data indicate no asymmetry. Ce, cerebellum; 
OB, olfactory bulb; P, pallium; poc, postoptic commissure; Pr, pretectum; SP, subpallium. 



unknown mechanisms generate subtypes among Otp- specified DA 
neurons, or different signalling environments of the anatomically 
distinct Otp groups direct different projection behaviours. 

Ascending DA projections. We analysed ascending DA projec- 
tions from the ventral diencephalon to the telencephalon. Only 
the more rostral Otp -dependent DC2 and DC4 neurons develop 
such projections. Among the 111 DC2/4 DA neurons analysed, ten 
somata (9%) sent an ascending projection to the SP (Figs 2, 6 and 7; 
Supplementary Fig. S2). Given that there are on average only 14.3 
(±2.02) DC2/4 DA neurons in a 4-d.p.f. brain, these data reveal that 



there are only one or two DC2/4 neurons with ascending projections 
per larvae (Supplementary Table S2). In addition, these large, pear- 
shaped unipolar cells 4 have fine local bifurcations and send proc- 
esses to hypothalamic domains, as well as far- descending axons 
that target the hindbrain and may traverse the spinal cord (Figs 6 
and 7; Supplementary Movie 4). Such projection characteristics are 
also typically observed for the majority of DC2, 4 and 5 neurons that 
do not send ascending projections (Fig. 7; Supplementary Fig. S6). 
Our counts of DC2/4 somata (left: 7.5 ±1.6, right: 6.8 ±0.9; 
Supplementary Table S2) are similar to those reported for the adult 
zebrafish brain 22 . This confirms that the number of large, pear-shaped 
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Figure 4 | Dopaminergic groups with characteristic short-range or local projection behaviours, (a-e) Z projections of confocal stacks in dorsal (Dor) 
and lateral (Lat) views, (a) A pretectal dopaminergic (DA) neuron (arrowheads) forms a local arbour within the pretectum (Pr) and projects into the 
ipsilateral tectum (T; magnified in right panel) and across the posterior commissure (pc) to the contralateral tectum (left panel). H, hypothalamus, (b) 
A diencephalic cluster (DC) 1 neuron (arrowheads) sends short proximal processes within the ventral thalamus (THv), into the endohypothalamic tract 
(eht) (magnified in middle and rightmost panels), and descending projections into the rhombencephalon (r). (c) A preoptic DA neuron (arrowhead) 
projects into the postoptic commissure (poc) and eht (magnified in right panel), (d) A DC3 neuron (arrowhead) located in the paraventricular organ 
(PVO) is liquor contacting (arrows pointing to end-foot) and contributes to the eht and hypothalamus (magnified in right panel; inset shows GFP signal), 
(e) Liquor-contacting DC7 soma (arrowhead) within the posterior recess possesses sparse processes (magnified in right panel) that target locally in the 
caudal hypothalamus (He). Anterior at left, dorsal at top; 4d.p.f. larvae; anti-TH (red), anti-GFP (green, white in panels with black contrast). Scale bars, 
20 jim. (f) Scheme of DA groups with local and intermediate projections (4d.p.f.). Projection pathways (lines) and target areas (arrowheads) are depicted 
in lateral (left) and dorsal (right) views. Short arrows and 1' denote local arbours or processes. An asterisk (*) marks liquor-contacting DA groups. For 
visual clarity, lines and arrows are not depicted left-right reciprocally in dorsal views. Our data indicate no asymmetry, ac, anterior commissure; AP, area 
postrema; LC, locus coeruleus; and MO, medulla oblongata. 



DC2/4 DA neurons does not increase from larval to adult stages. 
Analysing laterality, we found no preferential DA ascending pro- 
jections from either side of the posterior tuberculum (one-tailed 
paired f-test > 0.278, Supplementary Table S3). 

Prominent local projections of subpallial DA neurons. As the 

majority of telencephalic DA fibres do not originate from the ven- 
tral diencephalon, we examined the origin of the dense THir fibres 
within the SP. Among 60 GFP-labelled subpallial DA neurons ana- 
lysed, 55 (91.7%) were found to radiate dense local arbours within 
their territory. The prominent arbourization appears to cover the 
whole SP, whereas DA axons appeared not to project deep into the 
pallium. No projection to olfactory telencephalic territories was 
detected. Of the 60 somata, 25 sent a single contralateral process 
into the region of the anterior commissure. A substantial fraction 
of subpallial DA neurons also project a descending axon into the 
dorsal thalamic (12 of 60) or hypothalamic (18 of 60) regions, often 
in combination with local subpallial arbours and contralateral projec- 
tions (Figs 4f and 8; Supplementary Fig. S8, Supplementary Movie 5). 
Given the late differentiation of subpallial DA neurons, the dorsal 
thalamus may also be an intermediate target at 4 d.p.f. of hypothala- 
mus-bound descending projections. To determine whether left- right 
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asymmetry exists in the diencephalon-bound descending axons, the 
number of subpallial somata and projections on both sides of the tel- 
encephalon was counted. There are an average of 10.7 ( ± 1.9) and 10.3 
( ± 2.5) neurons on the left and right SP, respectively (two -tailed paired 
t-test > 0.379, Supplementary Table S4). We observed 14 descending 
processes on the left and 16 on the right SP (one-tailed paired 
f-test > 0.352, Supplementary Table S5). We therefore conclude that 
no chirality in the number of telencephalic descending connections 
to the diencephalon was observed. Semiquantitative abstraction 
of our tally of descending projections from subpallial DA neurons 
(30 of 60) predicts that in a 4 d.p.f. larval brain containing 21 ( ±3.9) 
of these somata, 10-11 neurons (50%) are likely to project to the 
diencephalon, therefore forming the major DA connections between 
the rostral and caudal forebrain. Thus, in quantitative terms, the pre- 
dominant DA connection between the SP and ventral diencephalon 
is not ascending, but descending. 

Discussion 

Understanding the projection behaviour of CA neurons is essential 
to evaluate their modulatory input into circuits that control behav- 
iours ranging from motor initiation to addiction. We used genetically 
mosaic zebrafish to systematically label each type of CA neuron, and 
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Figure 5 | Projection behaviours of otp-dependent dopaminergic diencephalic clusters, (a-d) Z projections of confocal stacks in dorsal (Dor) and lateral 
(Lat) views. Lateral views correspond to magnified regions (middle panels), (a) A DC2 neuron (arrowheads) contributes to the endohypothalamic tract 
(eht), forms local arbours and ventral processes into preoptic region (PO) and hypothalamus (H) (magnified in middle panel). A long axon descends 
into the hindbrain with branches (arrows) across the midline, and into the spinal cord (not shown), (b) A DC4 soma (arrowheads) with short processes 
contributing to the eht and connecting within the hypothalamus (magnified in middle panel), a dorsal branch extending towards the lateral pretectum 
(Pr) and tectum (T) (bottom panel), and long descending axons (arrows) into left and right hindbrain. (c) A DC5 soma (arrowheads) with ipsilateral and 
contralateral contributions to the eht (magnified in middle panel). From the ipsilateral eht, a branch (arrow) extends dorsoposteriorly into the medial 
longitudinal catecholaminergic tract to the hindbrain. (d) A DC6 axon (arrowheads) arbourizes extensively in the hypothalamus (magnified in middle 
panel) and projects towards the caudal midbrain and anterior hindbrain (bottom panel), r, rhombencephalon. Anterior at left, dorsal at top; 4d.p.f. larvae; 
anti-TH (red), anti-GFP (green). Scale bars, 20jim. (e) Schematic overview of DC2, DC4, DC5, and DC6 projection pathways (lines) and target areas 
(arrowheads) in lateral (left) and dorsal (right) views. 1' denotes local arbours or processes. For visual clarity, lines and arrows are not depicted left-right 
reciprocally in dorsal views. Our data indicate no asymmetry. AP, area postrema; dl, lateral diencephalon; LC, locus coeruleus; MO, medulla oblongata; 
PT, posterior tuberculum; SP, subpallium. 



to track their complete projections in three-dimensional (3D) space. 
This work generated the first complete view of the projection pat- 
terns that originate from CA somata in a larval vertebrate, provid- 
ing a close representation of the CA neurotransmitter class-specific 
projectome' 21 . There are two inherent limitations in this study: first, 
our focus on the larval stage at 4d.p.f. may not reveal aspects of 
the CA systems that mature later, including potential additional 
axonal tracts. However, important aspects of the larval DA systems 
resemble those found in the adult brain: 4 d.p.f. larvae have a diverse 



behavioural repertoire 3 , and the number of neurons in some CA 
clusters is highly similar (LC NA neurons 17 ) or is in the same order 
of magnitude (Otp-dependent DC2 and DC4 DA neurons 9,20,22 ) as 
in adults. Second, projections of £/z2-expressing neurons were not 
analysed, as our BAC and anti-TH antibody are specific for th and 
TH. However, several domains of th2 expression overlap with those 
of THir neurons in the diencephalon 9,10 , and most th2 DA neurons 
are caudal hypothalamic, which may likely possess mainly local 
projections similar to DC7 in zebrafish and A12 in mammals. 

7 
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Figure 6 | Individual somata of DC2 and DC4 emanate both ascending 
and descending projections, (a, b) Z projections of confocal stacks in 
dorsal (Dor) and lateral (Lat) views, (a) A DC2 neuron (arrowheads) 
within the posterior tuberculum (PT) projects both an ascending axon 
(single arrows) to the subpallium (SP) and a descending axon into the 
spinal cord. This neuron also forms a local arbour in the PT and ventral 
processes targeting the preoptic region (PO). Branches derived from 
the main descending axon were observed in the rhombencephalon (r), 
however, this region (double arrows) could not be clearly segmented 
because of the presence of a second GFP-labelled DC2 cell in the left 
PT, which contributed a midline-crossing branch into this area (inset). 
See Supplementary Movie 4. (b) A DC4 neuron (arrowheads) in the PT 
sends long ascending (single arrows) as well as descending (broad arrow) 
projections. Local processes in the PT and hypothalamus (H) also originate 
from this DC4 soma (double arrows). Its ascending projection arbourizes 
within the SP (magnified in the right panels). Anterior at left, dorsal at 
top; 4d.p.f. larvae; anti-TH (red), anti-GFP (green). Scale bars, 20jim. 
(c) Schematic overview of DC2 and DC4 projection pathways (lines) and 
target areas (arrowheads). Dorsal (top) and lateral (bottom) views. 'L' 
denotes local arbours or processes. DC2 and DC4 neurons contribute to 
wide-ranging ipsilateral and contralateral projections into all major brain 
subdivisions. For visual clarity, lines and arrows are not depicted left-right 
reciprocally in dorsal views. Our data indicate no asymmetry. AP, area 
postrema; dl, lateral diencephalon; LC, locus coeruleus; MO, medulla 
oblongata; Pr, pretectum; T, tectum. 
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Figure 7 | Diversity of projection behaviour subtypes of individual 
DC2 and DC4 neurons. DA neurons of one anatomical group express 
high diversity of projection behaviours within defined repertoires. DC2 
(69 somata) and DC4 (42 somata) neurons were analysed and shown 
to connect to tel-, di-, mes- and rhombencephalon. Individual somata 
usually connect through the endohypothalamic tract and send different 
combinations of projections to the subpallium, lateral forebrain, lateral 
tectum, hindbrain and/or spinal cord. DA innervation in the preoptic region 
and lateral forebrain primarily derives from DC2 cells. Individual somata 
were classified on the basis of their major projections (green), but smaller 
deviations in branching pattern and targets were allowed (yellow and light 
green). The rightmost column indicates the number of cells observed for 
each specific pattern (rows), and in parentheses the total number of cells 
analysed for DC2 and DC4. Dor, dorsal projections; lat, lateral projections; 
single in endohypothalamic tract indicates contribution of single axon to 
tract; (#) ipsilateral distal branch that crosses the midline to target the 
contralateral side. 



Our analyses of projection paths and target areas identified projec- 
tion behaviours typical for each CA neuronal group. However, within 
the repertoire of projections characteristic for each group, individual 
neurons exhibited high variability in their projection paths and tar- 
gets. The observed diversity may represent true additional genetically 
defined subtypes. Alternatively, projections formed by neurons of one 
subtype may depend on local guidance cues and may have stochastic 
aspects. Diversity in wiring diagrams of a single neuron type has pre- 
viously been shown for vertebrate motor neuron connectomes 23 , in 
contrast to the more stereotyped wiring in invertebrates 24 . 

Our study sheds light on evolutionarily conserved and dynamic 
aspects of CA systems in vertebrates. Compared with mammals 25 , 
our data revealed that all significant aspects of zebrafish NA pro- 
jections from the LC and MO appear to be conserved. Although 
projections from caudal rhomb encephalic NA groups into the fore- 
brain have not been reported for adult cypriniform fish 11,26,27 , other 
far-projecting brainstem transmitter systems have been described 
in perciform fish 28 . Our findings are consistent with a fundamental 
modulatory role of LC NA connectivity in global behavioural pat- 
terns, ranging from basic motor behaviours and stress responses to 
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mood and cognition 29 . The highly conserved nature of the LC NA 
system and its postulated integrative role in setting basic states of 
the brain may inform an understanding of neurodegenerative 
diseases of the CA systems, including Parkinsons disease 30 . 

The detailed analysis of DA ascending and descending systems 
in the forebrain surprisingly revealed that prominent DA connec- 
tivity between the telencephalon and diencephalon is contributed 
both by descending projections of subpallial DA neurons and by 
ascending projections from the ventral diencephalon. The posterior 
tuberculum was previously shown to be a source of ascending DA 
projections into the SP in zebrafish 11,12 . We extend these findings 



Figure 8 | Extensive arbourization of subpallial DA neurons in the 
telencephalon, (a-c) Z projections of confocal stacks in dorsal (Dor) and 
lateral (Lat) views, (a) Three GFP-tagged subpallial neurons (arrowheads) 
reveal extensive arbourization within the subpallium. Subpallial neurons 
also frequently connect the left and right subpallium through the 
anterior commissure (ac; magnified in right panel). Furthermore, one 
of the subpallial DA neurons sends a descending projection into the 
hypothalamus (H; arrows), (b) Two subpallial DA neurons (arrowheads) 
form dense local arbours but no distant projections (magnified in middle 
and right panels), (c) A single subpallial soma (arrowheads) radiates 
dense local processes ventrally into telencephalic territories (magnified in 
middle and right panels), projects in the anterior commissure across the 
midline to the right subpallium (SP), and sends a descending axon into 
the hypothalamus (arrows). See Supplementary Movie 5. Anterior at left, 
dorsal at top; 4d.p.f. larvae; anti-TH (red), anti-GFP (green; white in middle 
panels with black contrast). Scale bars, 20 um. 



by demonstrating that, at least by 4d.p.f., Otp-dependent poste- 
rior tubercular DC2 and DC4 DA neurons are the only source of 
ascending DA activity. These ascending connections are surprisingly 
scarce in the larval brain, with on an average one or two ascending 
connections on each side of the brain. Individual Otp-dependent 
DA neurons with ascending DA projections also connect with other 
major axon branches to diencephalic, hindbrain and spinal cord 
targets. Our observations suggest that Otp-dependent DA neurons 
may potentially be the only DA system that integrates ascending 
and descending circuits contributing to motor control. They possi- 
bly serve a role in neuromodulatory integration of sensory function 
(pretectum, tectum, thalamus), neuroendocrine state and physiol- 
ogy (hypothalamus), motor control (hindbrain and spinal levels), 
as well as control of behaviour and cognition (SP). Mammalian 
Otp-specified All DA neurons 15 have also been demonstrated to 
develop ascending neocortical and descending diencephalospinal 
projections 3133 , indicating that single-cell integration of rostral and 
caudal projections may be a component of an evolutionarily ancient 
DA neuromodulatory system. However, these studies showed 
neocortical/pallial but little or no subpallial All projections, raising 
the possibility that the final projection targets may not have been 
reached in 4-d.p.f. zebrafish larvae. Little is known about All func- 
tion, although All-lesioned rats are hyperactive 34 , and the system 
has been linked to Restless Legs Syndrome 35 and descending com- 
mand systems 36 . 

Given that All DA neurons also express the hypothalamic 
marker Nkx2.1 (ref. 15), an early evolutionary origin of this system 
is possible. An anterior DA cell population that sends projections 
ventrally as well as descending into the spinal cord has been shown 
in Amphioxus, although a homology assignment to a brain region 
in vertebrates has not been possible 37 . In the protochordate Ciona 
intestinalis, dopamine cells localize to a hypothalamus -related 
domain of the sensory vesicle 38 . Diencephalospinal DA projections 
have also been described in agnathan vertebrates (lampreys) 39 . All- 
type DA neurons from fish to mammals are cospecified with neu- 
rosecretory hypothalamic neurons 15,40,41 , which also develop dien- 
cephalospinal projections. This has led to the hypothesis that the 
vasopressin/ oxytocin systems may be involved in the spinal cord 
level control of sensorimotor, somatomotor and autonomic sys- 
tems 42,43 . We therefore hypothesize that the Otp-dependent DA and 
neurosecretory systems may be the ancestral systems that synergis- 
tically set basic behavioural patterns such as active (fight/flight) or 
passive (freezing) avoidance. 

Our data indicate that the previously proposed homology between 
the zebrafish posterior tubercular DA clusters and the mesodien- 
cephalic DA groups A8-A10 residing in mammalian substantia 
nigra (SN) and ventral tegmental area 44 46 needs to be reconsidered. 
Instead, the posterior tubercular DA groups should be considered 
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All homologues. This homology is not only based on similarities 
in projection pattern, but also on conservation of transcription fac- 
tors specifying these neurons 15,41 . Disputes about the transcription 
factor Nr4a2/Nurrl controlling terminal differentiation in mice SN 
and ventral tegmental area DA neurons but not zebrafish posterior 
tubercular DA clusters 44,47, further strengthen the notion that these 
clusters are different. 

The numerous dense arbours of the subpallial DA neurons 
we observed within the SP indicate that a predominant telen- 
cephalic dopamine source in zebrafish may be derived locally. In 
all vertebrates, the striatum is one of the subdivisions of the SP 48 . 
In zebrafish, the striatum is a dorsal subdivision of the dorsal SP 49 . 
However, because of the eversion process of the telencephalon in 
teleostei, and the lack of specific molecular markers, the delimita- 
tion of the striatum within the zebrafish SP has not been identi- 
fied. In contrast to the nigrostriatal contribution, dopamine input 
from intrinsic striatal dopamine neurons has received little atten- 
tion in mammalian systems. Striatal DA neurons, confirmed as 
TH and DA transporter expressing small bipolar neurons 50 , have 
been described in both primates 51 and rats 52 . Morphologically, the 
extensive local subpallial DA arbours in zebrafish are reminiscent 
of cortical arbours derived from a single neuronal projection from 
the SN in the rodent midbrain 53 . In contrast, the ascending projec- 
tions from o^p-dependent posterior tubercular DC2/4 DA neurons 
do not arbourize in the larval SP, and thus do not broadly connect. 
Therefore, it is conceivable that the predominant DA input into the 
SP in zebrafish is derived from intrinsic DA neurons. Our findings 
that subsets of zebrafish subpallial DA neurons emanate contralat- 
eral subpallial projections and send descending projections into the 
caudal hypothalamus, indicate that these neurons may contribute to 
circuitries modulating neural functions in distant brain territories. 
Cell-based regenerative therapies have so far focused on midbrain 
A9-A10-type DA neurons or precursors, to be ectopically placed 
into the striatum 54 . Our description of an endogenous subpallial 
DA system should encourage research into sources for endogenous 
striatal DA neurons, which may potentially be more stable in the 
striatum than SN- derived cells. 

In contrast to the highly conserved hindbrain NA systems, major 
aspects of the far-projecting DA systems appear highly dynamic 
during evolution. Given overlapping target area profiles, we specu- 
late that the NA and DA systems may have partially overlapping 
functions, as caudal hindbrain CA groups were reported to contain 
both NA and DA neurons in zebrafish and non-mammalian verte- 
brates 1,6 . Although the NA systems might have remained stable dur- 
ing the teleost to tetrapod to mammalian transitions, the evolution 
of the distantly projecting DA systems may reflect the development 
of higher order control through the sequential addition of neuro- 
modulatory mechanisms 55 . At an early stage, the local subpallial 
contribution may have been the predominant DA input into higher 
brain function. The increasing prominence of mesodiencephalic 
ascending DA projections may have facilitated integration of addi- 
tional inputs that are, for example, derived from sensory systems. 
This includes the progressive involvement of the cortex/pallium in 
processing thalamic sensory information relayed through the basal 
ganglia of tetrapods 56 . Whether the mesodiencephalic ascending 
DA systems evolved by caudal expansion of ventral diencephalic 
DA groups 56 or whether they evolved de novo in the context of 
mid-hindbrain organizer activity still needs to be resolved 57 . 

Methods 

Recombinant th:rasGFP BAC clones. We isolated a BAC clone (CHOR- 
B736O0777Q2) from the CHORI-211 zebrafish BAC library (http://bacpac.chon. 
org/zebransh211.htm). Using primers pTARBAC2.1_F (5'-TTGCGGCCGCTAAT 
ACGACTCACTA-3') and pTARBAC2.1_lR (5'-ATTGAGGCACTTGGTATTTT- 
TCGTT-3') for the ends of the BAC vector backbone, we sequenced both ends of 
the th insert to reveal that it contains 154kbp upstream and 12kbp transcribed 
region of the zebrafish tyrosine hydroxylase (th) gene (Supplementary Fig. SI). Two 



strategies were used to integrate ras membrane-anchored EGFP (mGFP) into the 
first exon of the th BAC clone. th:rasGFPsv (shuttle vector) was engineered using 
the R6Ky RecA-based shuttle vector pLD53.SCA-E-B modified to contain a 1-kbp 
fragment 5' of the th ATG and the mGFP reporter to be in-frame with the TH 
protein. Homologous integration into the BAC was verified by PCR (Supplemen- 
tary Fig. SI). thirasGFPsv was used in this study without removal of the vector 
sequence. 

thirasGFPvf (vector- free) was generated to reduce potential effects due to 
vector sequences, by using galK selection based recombineering 58 . Amplification 
of the galK selectable marker from pGalK was carried out using forward primer 
5'- CCAGCGCTCCCAAGAGACCACAATATAGTTCCATCTGTTCCGCGCGCG 
CGCCATC CCTGTTGACAATTAATCATCGGCA-3' and reverse primer S'- GTAA 
TAGAATCTGAGCGCTCGAGTTCAGACGCAGCTCTCCGGATGCTCTTTGT TC 
AGCACTGTCCTGCTCCTT-3' (th homology arms are in italics and underlined). 
The PCR product was electrop orated into transformed bacteria and recombinants 
selected from Gal indicator plates were used in the second recombineering step. 
Forward primer 5'- CCAGCGCTCCCAAGAGACCACAATATAGTTCCATCTGTTC 
CGCGCGCGCGCCATC CCACCATGGTGAGCAAGGGC-3' and reverse primer 
5'- GTAATAGAATCTGAGCGCTCGAGTTCAGACGCAGCTCTCCGGATGCTC 
X7TGIAAAAAACCTCCCACACCTCCCC-3' were used to amplify the targeting 
cassette from mGFP (sequences are in bold) using proofreading Taq Polymerase 
(Advantage Taq, BD Clontech). Recombinants from Gal counter-selection plates 
were confirmed by PCR and sequencing. thirasGFPvf was subsequently trans- 
formed from the SW102 Escherichia coli strain back into DH10B host cells. Double 
acetate precipitation and caesium chloride gradient separation were carried out 
to purify the modified BACs. DNA quality and concentration were checked by 
pulsed-field electrophoresis. Restriction digests, Southern hybridization, sequenc- 
ing and PCR characterizations were carried out to verify single copy recombination 
of the reporter gene in the correct orientation (Supplementary Fig. SI). Compar- 
ing results of CA neuron labelling in embryos for both BAC constructs, we have 
observed no undesirable backbone effects for thirasGFPsv, which contributed to 
81.6% of our analyses. 

Random GFP-labelling of tfi-expressing neurons. Zebrafish were main- 
tained and bred under standard conditions at 28.5 °C. For microinjection into 
the cytoplasm of one-cell stage embryos, the modified BACs were prepared at 
100-120ngjil _1 in 0.1 M potassium chloride and 0.05% phenol red. To avoid mela- 
nin pigment formation, embryos were incubated in 0.2 mM l-Phenyl-2 thiourea 
(Sigma). At 3-4d.p.f., embryos were mounted in 2.5% methylcellulose containing 
Tricaine and screened for individual GFP-labelled neurons. 

Immunolabelling and confocal microscopy. Immunofluorescence and immuno- 
histochemistry were carried out as described 12 using rabbit anti-tyrosine hydroxy- 
lase antibody (1:500 (ref. 15)) and chicken anti-GFP (1:400, Molecular Probes) or 
mouse anti-GFP antibodies (1:1,000, Roche). Alexa555- or Alexa594-conjugated 
goat anti-rabbit IgG (1:2,000 and 1:1,000, respectively, Molecular Probes) and 
Alexa488-conjugated goat anti-chicken IgG or anti-mouse IgG (1:600 and 1:1,000 
respectively, Molecular Probes) were used as secondary antibodies. Embryos were 
mounted with 80% glycerol/ 1.2% low-melting agarose/PBS. 

Confocal Z stacks were recorded with a Zeiss LSM 510 laser- scanning confocal 
microscope with pinhole diameters set to one airy unit for Alexa555 and Alexa488- 
labelled embryos using the 488- and 561-nm laser lines. Embryos stained with 
Alexa594 and Alexa488 were examined under Zeiss LSM 510 NLO with pinhole 
diameters between 360-520 Jim using the two-photon 770-nm laser line. All stacks 
were recorded with a 25x NA 0.8 multi-immersion lens (LD-LCI Plan-Apochro- 
mat) and glycerol immersion medium. Stack sizes are 1,024x1,024 pixel in xy 
direction, with the number of z planes depending on the thickness of the sample. 
The voxel sizes used were 0.5 umx0.5 Unix 1.3 jam in the z direction. Pixel times of 
1.26-3.20 (is were used. 

To better visualize neuronal projections in the context of the anatomical struc- 
tures, differential interference contrast microscopy was used (Supplementary Fig. 
S9). Following confocal analysis, embryos of interest were incubated overnight in 
1% goat serum/PBST containing DSB-X biotin-conjugated goat anti-chicken IgG 
(1:1,000, Molecular Probes), processed using biotinylated HRP and avidin reagents 
(Vectastain ABC HRP kit), and stained with the chromogen diaminobenzidine. 
Images were documented with DIC Axioscope2/AxioCam MRc (Zeiss) and proc- 
essed with ImageJ (http://rsbweb.nih.gov/ij/). 

Image processing and analysis. For far-projecting neurons, adjacent whole-brain 
image stacks were recorded. To enhance the visualization of GFP-tagged somata 
and their projections, we sought a simple solution to isolate them from nonspecific 
signals. In the whole-larvae preparations, immunofluorescence noise, auto-fluores- 
cence (for example, some pigment cells), as well as other GFP expressing cells build 
up significant nonspecific signals when visualizations of image stacks are prepared. 
We wanted to avoid established noise-reduction approaches including thresholding 
or despeckling, given that fine CA processes in our stacks often had cross-section 
areas of very few pixels only. Therefore, we developed an image-processing pipeline 
for interactive in silico feature extraction, which is amenable to automated batch 
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operation and includes a quality control element that allows independent and 
objective verification of 3D objects (Fig. 1). 

The data set for a single brain recording was processed as follows (Fig. lb): 
Maximum intensity projections (MIP) in dorsal and lateral directions and MIP 
movies of the revolving 3D brain image were created for quick visualization. The 
user proceeded to draw the dorsal and lateral segmentation masks over the MIPs 
using our ImageJ HDF5 plugin (http://lmb.informatik.uni-freiburg.de/lmbsoft/ 
imagej_plugins/). The 3D intersection M of the 2D lateral L and dorsal D masks 
was computed as M(x,y,z) = D(x,y) L(x,z), where M, D and L are scalar fields 
containing 1 for object pixels and 0 for background pixels. The resulting 3D mask 
was applied to the green channel, thereby allowing the segmentation of labelled 
neurons within a single brain and disregarding nonspecific signals and background 
(Fig. lc-e). For quality control procedures as well as comparison to other projec- 
tion tracking approaches, see Supplementary Methods and Fig. S10. All figures 
were assembled using Adobe Photoshop CS2 or Adobe Illustrator CS4. Brightness 
and contrast levels were adjusted in Photoshop to ensure that fine axonal and 
dendritic processes appear in print. 

Identification and tracking of CA projections. Three-dimensional renderings 
of whole-brain image stacks were carefully analysed for the collection of complete 
information on individual neuron soma and projections. 

Step 1: Identification of labelled CA neurons. On the basis of their anatomi- 
cal position and correlation with anti-TH stained CA clusters, the neuronal CA 
subtype of the GFP- labelled CA neurons was identified. For each neuronal soma 
to be included in further analyses, the anti-TH (red) and anti-GFP (green) signals 
were analysed at high resolution and both colour channels were adjusted to verify 
the CA phenotype. For each neuron documented in the figures of this paper, the 
corresponding somata are displayed in Supplementary Fig. SI 1. 

Step 2: Evaluation and classification of CA projections. Using the location of 
THir somata and axon tracts as anatomical reference, the projections of individual 
labelled neurons and their target areas were identified and classified. Histologically 
stained experimental whole-mount larvae preparations also enabled us to correlate 
projection paths and target areas with anatomical information (examples in 
Supplementary Figs S3 and S9). 

Selectivity of th:rasGFP BAC labelled neurons. We assessed whether tkrasGFP 
BAC may drive ectopic expression of GFP in cells other than CA neurons. The 
embryos (9%) analysed by confocal microscopy also showed anti-GFPir cells that 
did not correlate with anti-THir (data not shown). However, these cells appeared 
randomly distributed, including cells outside the CNS, such as somite cells. Such 
events may be explained by stochastic breaks of the BAC during integration and 
enhancer trap effects. The ectopic expression did not interfere with our analysis 
of CA neurons, which in each case could be clearly identified based on anti-THir. 

Statistical analysis. To evaluate the likelihood that we may have missed signifi- 
cant projection patterns, we carried out statistical analysis. Given that different 
projection patterns for a specific neuron group may not be randomly formed and 
labelled, we analysed graphically how many projection patterns or target areas with 
a specific case number were identified for each DA group (Supplementary Tables 
S6, S7). Statistical analysis was carried out under the assumption that each pattern 
occurred at random, which is a simplification given that different groups had dif- 
ferent cell numbers and differences in complexity of projection patterns. Binominal 
regression analysis was used to identify the interception with the y axis. Plotting 
the number of times a specific projection pattern or target area was observed, 
revealed that we may have missed about 10% of potential target areas and about 
20% of projection types. However, the several-fold increase in number of analysed 
embryos required to reduce the zero class significantly was not technically feasible 
in this study. 

A non-random effect on the level of saturation would hypothetically be possi- 
ble if the th:rasGFP BAC has a bias in labelling only a subset of DA neurons. How- 
ever, we observed that each group contributes to the total number of GFP-labelled 
cells at approximately the same ratio as its THir cell count compared with the total 
number of THir cells. Therefore, we believe that we have surveyed the complete set 
of THir cell types in the larvae. However, we cannot fully exclude the underrepre- 
sentation of a specific THir cell type. 
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